Tunable adsorbate-adsorbate interactions on graphene 
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We propose mechanisms to control the interaction between adsorbates on the surface of graphene. 
We show that the interaction between a pair of adsorbates — the change in adsorption energy of one 
adsorbate in the presence of another — is dominated by a non-linear exchange interaction mediated 
by graphene's 7r-electrons. A b initio density functional, numerical tight-binding, and analytical 
calculations are used to develop the theory. We demonstrate that the interaction can be tuned in a 
wide range (up to complete reversal) of values as a function of distance by adjusting the adsorbate- 
graphene bonding or the chemical potential. 



Graphene is a one- atom- thick sheet of carbon. It is 
of interest for a variety of applications ranging from 
transport to photovoltaics and DNA manipulations []]- 
and has potential to become a successor to silicon in 
micro- and nanoelectronic applications. The unique elec- 
tronic transport properties of graphene are rooted in its 
two-dimensional honeycomb carbon network structure. 
However, wider applications are limited by difficulties in 
opening a bandgap and its lack of processability. Both 
of these limitations can be addressed through chemical 
modifications Recently, the range of potential ap- 
plications has been enlarged by recognizing that chem- 
ical functionalization [5|4Zl of graphene can change its 
electronic, chemical, and mechanical properties, tailoring 
them to meet specific needs. As a result, design and con- 
trol of functionalization has become an important chal- 
lenge. 

Atoms or molecules near the surface of graphene can 
adsorb via Van-der-Waals interaction (physisorption) or 
by forming chemical bonds with one or several car- 
bon atoms on the surface (chemisorption). Chemi- 
cal functionalization of graphene is typically based on 
chemisorbed atoms or molecules that remain chemically 
active or provide other functions, e.g., influence con- 
ductance when adsorbed in large numbers. Adsorbate- 
graphene interaction for variou s ty pes of adsorbates have 
been discussed in literature [8l-ll6|. Recently, it has been 
recognized that adsorption energy of an adsorbed atom 
or molecule can depend strongly on the presence of other 
adsorbates in its vicinity. The adsorbate-adsorbate (A- 
A) interaction becomes especially important at large con- 
centrations of adsorbates, affecting both stability and ag- 
gregation when forming functionalized surfaces [17l-l2lj. 
Therefore, understanding the nature of the A- A interac- 
tion is a key to design and control of functionalization of 
graphene. 

In this Letter we investigate the microscopic mecha- 
nisms underlying the interaction between two atoms or 
molecules adsorbed on graphene at a distance from each 
other. We demonstrate that A- A interaction is domi- 
nated by an exchange interaction mediated by graphene's 
7r-electrons. We show that A-A interactions can be modi- 



fied up to complete reversal of sign by changing the chem- 
ical potential or the adsorbate-graphene bonding. We use 
a combination of first-principles density functional theory 
(DFT) calculations [22|, numerical tight-binding (TB), 
and analytical functional integral approach to identify 
the role of different microscopic mechanisms responsible 
for the interaction. We illustrate the approach using flu- 
orine adsorbates as an example relevant to a variety of 
applications, and demonstrate how our analysis applies 
to more complex systems. 

Chemisorbed atoms or molecules can form one (mono- 
valent) or several (multivalent) bonds with the carbon 
(C) atoms of the graphene lattice. Monovalent adsor- 
bates, e.g., fluorine (F), adsorb on top of one of the C 
atoms forming (typically) a mixed covalent-ionic bond. 
In the case of F the covalent component of the bond is 
formed by the overlap of F pz-orbital (z-axis perpendicu- 
lar to the grahene's plane) and s-pz hybridized orbital of 
C atom. The C atom participating in the bond changes 
its hybridization from sp 2 to sp 3 taking one pz electron 
out of the 7r-system. This does not completely remove 
the amplitude for the electron to hop to the fluorine- 
carbon complex, but instead it leads to creation of a more 
complex state near each adsobate site as well as a local 
deformation of graphene. These modifications can lead 
to four possible A-A interaction mechanisms: (i) direct 
overlap of adsorbate's electron orbitals if two adsorbed 
atoms are in direct contact with each other; (ii) gain or 
loss of energy due to change of elastic deformation energy 
of the graphene sheet when two adsorbates are nearby; 
(iii) Coulomb interaction between excess charges drawn 
to each adsorbate; (iv) interaction induced by multiple 
scattering of graphene's 7r-electrons off the adsorption 
sites. The first contribution is effective only for atoms 
adsorbed in direct proximity from each other [23]. Lat- 
tice deformations can be demonstrated to produces 1/r 3 
contribution (see also 24], HH) that is weak ~ 1 — 10 
meV, except at the shortest distances. In what follows 
we demonstrate that the primary contribution to A-A 
interaction at larger distances comes from the interplay 
between Coulomb and scattering-induced, or electron ex- 
change, interactions. 
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The Coulomb and electron exchange interactions in the 
system of two adsorbates on graphene can be described 
by a TB Hamiltonian [28[ of the form 

(ij)' n=l,2;s (in),s 
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Here the first term describes graphene's 7r-electrons 
26L 27 1, the cj (cj) are creation (annihilation) operators 
of states on C atoms [29[ , and t is the hopping inte- 
gral between the nearest pz states. The adsorbates (num- 
bered by n = 1,2) can form chemical bonds with several 
nearby carbon atoms removing their pz states from the 
graphene's 7r-system, as indicated by the prime in (ij) f . 
Each bonded adsorbate-carbon complex (ACC) can have 
several states, as indicated by s in the second term of 
(PP). The hopping integrals between ACC states and pz 
states on other C atoms is denoted by t' n s in the third 
term of (pQ). The creation (annihilation) operators for s- 
th ACC state are a) n s (a n>s ). The fourth term of (pQ) is the 
Coulomb interaction between the 7r-electrons (n = c^c). 
The fifth term represents the Coulomb interaction be- 
tween the 7r-electrons and electrons that occupy ACC 
states (ft = a^a). Finally, the sixth term is the Coulomb 
interaction energy between electrons on the two ACCs. 

The parameters u and t f of the Hamiltonian (pQ) are 
not known. They can be obtained from experiment or 
from DFT calculations [lH, l30l - [33| by analyzing the den- 
sity of states (DOS) of graphene with a single adsorbate. 
We will consider F as the adsorbate to illustrate this 
parametrization. Fluorine bonds to one carbon atom. 
The two states that are closest in energy to graphene's 
Dirac point and are sufficiently connected by tunneling 
to 7r-electrons are the bonding (symmetric) and the anti- 
bonding (antisymmetric) overlap of pz orbitals of C and 
F. These states can be identified in the projected DOS 
(PDOS), as shown in Fig. QJa). In the case of F, it is 
suffisient to consider only the anti-bonding state in the 
TB description. In this case s in Eq. (pQ) takes a single 
value an will be omitted. 

To extract the values of u and t' we compute the TB 
DOS [34] for one adsorbate based on the first three terms 
of (pp). The resulting DOS [see Fig. [TJb)] has two peaks 
associated with the adsorbate: one near the Dirac point 
and one at higher energies ~ u. The former is due 
to the "missing" (sp 3 -hybridized) carbon atom which is 
now a part of ACC. The non-zero tunneling between this 
state and the state of the C-F complex at ~ u shifts 
their energies away from one another. This shift can 
be obtained analytically in the limit of linear disper- 
sion relation as e(e — u) = 7r\/3£ /2 /2 log |3 1 / 4 v / 7rt/e|. For 
small t f , the two energy branches can be approximated 
as e « u ± yju 2 + (3t'/2) 2 /2 which can be understood on 
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Figure 1: The results of DFT, TB, and analytical calculations, 
(a) DFT DOS and projected DOS (PDOS) for the system of 
graphene with a single adsorbed F and clean graphene with 
one missing C atom. The bonding and anti-bonding states of 
the F-C complex can be identified as buildup of PDOS (pro- 
jected on pz states of F and C) above and below the Dirac 
point. The inset magnifies the Dirac point region and shows 
the shift of the DOS peak due to non-zero t' . (b) TB DOS as 
a function of t' . The two solid (blue) curves are the analytical 
solution to the energies of the adsorbate states as discussed in 
the text. The white curves are total DOS (the reference point 
for each curve is shown by a vertical line). The background 
density plot represents the difference between TB DOS with 
and without F atom: the two main peaks (darker shade) fol- 
low the analytical solution for small-to-intermediate values of 
t' It. (c) Interaction energy between two F's as a function of 
neighbor order distance as shown in panel (d), see the text 
for discussion, (d) Numbering convention for the distance 
(positions) between the two adsorbates: one adsorbate is at 
position "0" (red) and the other is at any of the other marked 
sites. 



the basis of a simple model of two states connected via 
tunneling amplitude t' . In order to track the shift due 
to t' we compute DFT DOS [30[ for graphene sheet with 
one missing carbon atom, which corresponds to u — )• oo 
and t' = 0. Comparing DFT and TB results we obtain 
t' ~ 0.27t and u ~ 0.75t. The parametrization for other 
types of adsorbates can be done similarly. 

The energy of the system of two adsorbates on 
graphene, £(R), described by (pQ) as a function of the 
distance R = R2 — Ri between the two adsorbates 
can be separated into three components: the energy of 
graphene with two adsorbates at R = 00, £ (00), the en- 
ergy of Coulomb interaction between adsorbates, E'c'(R), 
and the electron scattering-induced interaction energy 
E e (R), i.e. £(R) = £(00) + E(R), where the A-A in- 
teraction energy is E(R) = i?c(R) + ^e(R)- In the 
case of a simple adsorbate, such as F, the energy of 
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Figure 2: Adsorbate-adsorbate interaction energy for different 
scattering strengths and distances, (a) Crossover from strong 
to weak coupling limits with distance: at short distances the 
minima occur when adsorbates are on different sub-latices 
[odd sites along AC, in Fig. Hid)]; at R c of about AC=10 
(shaded area) the phase abruptly changes, and the minima 
occur when adsorbates are at the same sub-lattice sites (even 
site along AC). The dashed line is the weak coupling limit, 
(b) The A-A interaction as a function of t' and u at NO=3 
distance between the adsorbates. The parameters for two F 
adsorbates are indicated by an arrow and "F-F" . The bottom 
sub-panel shows t' = t dependence (solid curve), which start 
quadratically in u (dashed curve). The right sub-panel shows 
interaction energy for AC=5 and AC=11. (c) The A-A inter- 
action energy as a function of distance for different chemical 
potentials fi/t = 0,0.15,0.3,0.5: solid (red), dashed, dash- 
doted, and doted (blue) lines respectively. The top sub-panel 
shows the values of \i relative to graphene's spectrum. 



Coulomb interaction can be well approximated [35| by 
Ec(R) = e 2 Z 2 / Aite^r, where Z is the effective charge 
drawn to the adsorbate complex obtained from DFT cal- 
culations (Z « —0.39 for F). As we will see shortly, the 
screening due to the fourth and the fifth terms (pQ) do not 
alter the 1/R dependence in Ec(R>) and can be ignored. 
As an example, we consider the system of two F adsor- 
bates. To examine Ec(R) we performed a DFT calcula- 
tion [30] in which C atoms were restricted to their orig- 
inal (clean, fiat graphene) positions, while the positions 
of the F atoms were allowed to relax [see "DFT, flat" 
curve in Fig. (He)]. Each F attracts negative charge. The 
restriction (to flat graphene) enforces sp 2 hybridization 
for the C atom under each F atom making the C-F bond 
mostly ionic. It also eliminates possible elastic deforma- 
tion of graphene sheet around both F atoms, as well as 
electron scattering due to the second and the third terms 
in (pQ). The C pz orbitals under F's remains part of the 
graphene's 7r-system, i.e., t' n0 = t and u n $ —> 0. The 
electron exchange that can originate from the fifth term 
of (pQ) is still present in this configuration. However, by 



comparing Ec(R>) and "DFT, flat" results in Fig. []Jc) 
we conclude that it does not contribute substantially to 
the A-A interaction, even for the high values of Z in this 
example. 

We now investigate the role of the u and t' terms in 
E e (H). We compute the interaction energy E(H) be- 
tween to F's with fully relaxed DFT [30[, in which we al- 
low deformation of graphene ["DFT" curve in Fig. 0Jc)]. 
We also evaluate the TB [34| interaction energy E e (TV) 
based on the first three terms of (pQ). The TB E(R) is 
denoted as "TB with Coulomb" in Fig. DJc). The val- 
ues of u and t' found earlier for a single F were used. 
The DFT and TB curves show remarkable agreement at 
distances starting from the third nearest neigbour order, 
NO=3, [see Fig. 0Jd)]. Larger deviations at NO=l and 
NO=2 distances can be attributed to a combination of 
local deformation in graphene-adsorbate system and to 
direct overlap of fluorine electron orbitals not accounted 
for by Hamiltonian (pQ). To understand the physics of 
E e (TV) we will examine an analytical solution for this 
system in two limiting cases. 

The interaction energy E e (H) can be obtained by sum- 
ming all R-dependent diagrams for the electron free en- 
ergy using linked cluster expansion [36[. The expression 
for E e (R) can be rearranged to 

£ e (R)=2 J dYdY'Vi{Y)V 2 {v')J(R + Y-Y'), (2) 

J( r ij) = Jijfaj) = J ^Glfiiuj.Yi^Gjiiiuj.YjYi). (3) 

Here V\{2) is the scattering potential due to the second 
and third terms of Hamiltonian (pQ), r^- = — yj is 
the distance between two lattice sites, one at sub-lattice 
i and the other at j. The indexes i,j take the value 
of A and B when is on two different sub-lattices of 
graphene 0, E3. The function G ( °\ioj,r) is the bare 
equilibrium Green's function [37| and Gij(iuj,YiYj) is the 
dressed equilibrium Green's function which contains mul- 
tiple scattering events [36]. The limits u = 0, t' —> and 
t' = t, u —> oc are equivalent and represent a strong 
scattering system: in both cases an electron cannot en- 
ter the adsorbate site(s). A weak scattering limit can be 
illustrated by the system with u <C t, t' = t. We will 
focus on the t' = t case and set V n (r) = uVcS(y — R n ), 
where Vc is the volume of graphene's primitive cell. We 
are interested only in G(iu), R1R2), which can be found 
exactly 

G(zcj,RiR 2 ) = W(iu,$(iu,R))G°(iu,K)W(iu),0), (4) 

where Vl(iu,X) = [1 - uV c G°(ioj, 0)/y/2 - X}' 1 and 
£(zcj,R) = (uVc) 2 G (i(j,K)W(iu,0)G (iu,-Ii)/2. Ma- 
trix products are implied where r/'-indexes are dropped. 

When the scattering is small, u — >> 0, G — >> G^ and we 
obtain E e (R) = 2u 2 Vg J^ KKY (R), recovering the stan- 



dard RKKY range function for graphene [37J 

tRkky/tj \ J-iy- j+1 V3 l+co S [R(K-KQ + (2fl+7r)^] 
ij [ ' 3^2 6 7rfa 4 {R/af 

(5) 

where we use A — A = and A — B = 1 for z, j indexes. 

Unlike magnetic scattering in RKKY, which is weak 
in most cases, scattering of electrons on adsorbates can 
be strong. In the limit u — >• oo, the frequency inte- 
gral in Eq. (j3j) is given predominantly by w ~ v/i?. 
As a result, at large (but finite) distances we can use 
y\{iuj,${iLJ, R))9t(io;, 0) ~ — l/#(io;,R) and obtain 

(-3)*--7'7ta/3£ 1 + cosR(K-K') 
ij( ^ 2 6 (i?/a) 1+cos[R(K-K / ) + (2^+tt)%] ' 

(6) 

In some cases # ~ 1 + cos[R(K - K') + 2(9 + tt] 
and the frequency integration is determined by the renor- 
malization factor 9t(icj,0) 2 . In this case J^a(R) ~ 
— y/%ta\Y + cosR(K — K')]/4R. The interaction energy 
given by © and ([2]) is shown in the inset of Fig. []J along 
with with the numerical TB result for two F's. Devia- 
tions of the analytical result at short distances are due to 
approximate integration and the use of linear dispersion 
in the analytical estimates. 

The interplay between the two distinct regimes, (i) 
weak coupling from Eq. (j5j), and (ii) strong coupling from 
Eq. (|6]), can be used to modify the adsorbate-adsorbate 
interaction. Note that the A-A interaction in these lim- 
its have opposite phases, c.f. prefactors (— l)* - -^ 1 and 
(— 3)* _J . At finite scattering strength and R — >• oo the 
the range function J in Eq. (j3]) approaches the weak 
coupling limit. This can be most easily seen in the 
TB calculations for two F adsorbates placed at AC sites 
[see Fig. [2fa)]. The crossover between the two limiting 
regimes takes place gradually over a substantial distance. 
However, the phase associated with sub-lattice oscilla- 
tions noted above changes abruptly at some Rc [see the 
shaded region in Fig. Efa)]. The value of Rc can be es- 
timated from the denominator of 9t(za;, $(ioj, R)). We 
obtain Rc I a ~ u/t. 

To investigate the interplay of weak and strong cou- 
pling we perform TB calculations [34| for different values 
of t f and u. We use the sign of the A-A interaction at a 
specific distance R, e.g. NO=3, as an indication of the 
interaction regime [see Fig. Efb)]. The weak interaction 
regime is observed at small-to-intermediate u and for a 
range of t' values (depending on R). At large values of u, 
for t' — >• £, and finite distances the interaction converges 
to the strong coupling limit [see Eq. [6] as expected. The 
region with t' and u, which is typical of covalently bound 
adsorbates [indicated (approximately) by a shaded box in 
Fig. Efb)], is near the transition boundary. As a result, 
small changes in the bonding characteristics, e.g., due to 
changes in internal structure of an adsorbed molecule, 
can potentially lead to dramatic modification of the A-A 



interaction at different distances. This effect, combined 
with the recently proposed directed motion of the ad- 
sorbates [15] can provide control of interest in a variety 
of applications, such as chemical sensing and patterning. 
Ionically bound adsorbates are expected to leave the car- 
bon 7r-conjugated network intact as t' — >• t, interacting 
predominantly in the weak coupling regime. The change 
of sign between the strong and week coupling regime will 
also take place for multivalent adsorbates, e.g., oxygen 
(l5| . In that case the oscillatory behavior of the A-A 
interaction in the weak coupling limit will be similar to 
that for magnetic RKKY interaction when spin impuri- 
ties reside in between the nodes of the lattice [ItJ • In the 
strong scattering limit the renormalization factors 91 will 
change the range function as it happened in the exam- 
ple of the monovalent adsorbates above. As a result, the 
interplay between strong and weak coupling limits can 
be utilized to control the A-A interaction in a similar 
fashion. 

Finally, we investigate the role of the chemical poten- 
tial, \i. Changes due non-zero \i in 

G (0) G (0) 

are not signif- 
icant up to very high chemical potentials where the scat- 
tering processes are no longer confined to well-defined 
Dirac cones [see the top sub-panel in Fig. Efc)]. The 
renormalization factors in ((4]), however, are more sensi- 
tive. As a result, chemical potential can be used to alter 
the A-A interaction in the strong coupling regime. As 
an example, we consider the system of two F adsorbates 
as before. In Fig. [2fc) we plot E(H) for several values 
of ji from to 0.5t that can be achieved by back-gating 
or chemical dopping. Substantial changes in the A-A 
interaction can be seen even for relatively small values 
of \i\ e.g., in the example in Fig. [2jc) the minimum of 
the interaction energy shifts from NO=3 to NO=6 as \i 
increases. 

We have demonstrated that the dominant contribution 
to the A-A interaction is due to scattering of 7r-electrons 
on the distortions cased by the adsorbates. The proposed 
tight-binding model gives excellent agreement with the 
DFT calculations. Mechanisms to control the A-A inter- 
action at different distances has been proposed. 
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